A family of spin-S chain representations of SU(2)fc Wess— Zumino— Witten models 
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We investigate a family of spin-S 1 chain Hamiltonians recently introduced by one of us [1]. For 
5* = 1/2, it corresponds to the Haldane-Shastry model. For general spin S, we find indication that 
the low-energy theory of these spin chains is described by the SU(2)fc Wess-Zumino- Witten model 
with coupling k = 2S. In particular, we investigate the S = 1 model whose ground state is given by 
a Pfafnan for even number of sites N. We reconcile aspects of the spectrum of the Hamiltonian for 
arbitrary N with trial states obtained by Schwinger projection of two Haldane-Shastry chains. 
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Introduction — Conformal field theory (CFT) has sig- 
nificantly deepened our understanding of quantum spin 
chains and phase transitions in two-dimensional space- 
time [2 4]. In particular, it has provided a platform to 
distinguish different universality classes of critical spin 
chains according to e.g. the power law parameters of 
multi-point spin correlation functions. Moreover, the 
central charge parameter c of the CFT emerges in experi- 
mentally connected quantities such as specific heat [5, 6]. 
For critical SU(2) spin-S* chains, the conformal symme- 
try is supplemented by a Lie group symmetry and yields 
a low energy field theory description in terms of Wess- 
Zumino- Witten (WZW) models [7, 8]. The WZW action 
consists of a non-linear a term plus the Wess-Zumino 
term, whose topological coupling factor k is constrained 
to be integer and defines the level k of the WZWfe 
model [4]. In the absence of broken continuous symme- 
tries, WZWi is the generic low energy theory for half- 
odd-integer antiferromagnetic spin chains [9], while inte- 
ger spin chains show a Haldane gap [10]. Spin chains 
associated with WZWt>i models only appear at rare 
multi-critical points, and are of particular interest for in- 
teger spin chains where they can mark phase transitions 
between different gapped phases [9, 11, 12]. 

The 1/r 2 Heisenberg spin- 1/2 chain independently 
found by Haldane and Shastry [13, 14] plays a special 
role among all other models related to WZWi. While 
any finite size chain generally exhibits logarithmic cor- 
rections to the WZW long-wavelength limit depending 
on the system length L = aN where N is the number of 
sites and a is the lattice constant, the Haldane-Shastry 
model (HSM) exactly obeys the WZWi scalings for any 
finite L. (This is related to its enlarged Yangian quan- 
tum group symmetry [15].) Moreover, the exact HSM 
ground state wave function corresponds to the Laughlin 
wave function [16] in terms of bosonic spin flip particles 
at v = 1/2 filling, and establishes the notion of fractional 
quantization and statistics of spinons as the fundamental 
excitations of quantum spin chains [17]. In this respect, 



it is the one-dimensional analogue of the chiral spin liq- 
uid [18, 19], where the concept of topological order was 
introduced in spin models [20]. 

It was realized recently that not only the Laughlin 
state, but a subset of the bosonic Read-Rezayi quan- 
tum Hall series [21] can be generalized to a series of 
singlet spin S wave functions at spin flip particle fill- 
ing v = S [22]. For one spatial dimension, these states 
are constructed by a Schwinger boson projection tech- 
nique [23] of k copies of HSM ground states. The possi- 
bility of defining SU(2) invariant spin S = k/2 states 
related to the parafermionic CFT construction of the 
Read-Rezayi states is intuitive, as their current algebra 
reduces to the SU(2)fc Kac-Moody algebra [21] for the 
relevant fillings. As one particularly interesting mem- 
ber, the k = 2 state corresponds to an S = 1 Pfafnan 
spin state, which promises to establish the basis at which 
manifestations of non-Abelian spinons in spin chains can 
be investigated. 

From the combined view of state properties and low en- 
ergy theory, it is then natural to ask whether these states 
may establish finite size representations of spin chains as 
the conformally invariant fixed point of WZWfc>i in the 
same way as the HSM does for WZWi . A step in this di- 
rection has been recently accomplished by one of us [1], 
who introduced a family of Hamiltonians which singles 
out the spin S chain states obtained by projection from 
k = 25 HS ground states as exact ground states. In 
this Letter, we further investigate these Hamiltonians. 
We numerically show that these S = k/2 spin chains are 
critical and indeed connected to the WZWfe in the long 
wavelength limit. In contrast to the HSM, however, we 
find that the S > 1/2 models exhibit logarithmic correc- 
tions and hence do not describe the conformally invariant 
fixed point of WZWfc>i. For the S = 1 chain, we further 
analyze the excited states of the model, and generalize 
the Schwinger boson projection method to trial states 
for the simplest excitations. 

Hamiltonians and ground states — The general Hamil- 
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FIG. 1. (color online) Entanglement entropy S(x) of Hs 
for 5 = 1 and S = 3/2 vs. the conformal length d{x) = 
L/n sin nx/L. Data points are given for L — 20 (red), L = 24 
(blue), and L = 30 (green). The thermodynamic WZW re- 
sults are sketched by dashed lines. The fitted central charges 
(solid lines) agree well with SU(2) WZW2S field theory within 
reasonable error bars. 




FIG. 2. (color online) Scaling dimensions for the Haldane 
Shastry model (left) and Hi (right). Left: The finite size data 
(crosses) exactly matches the estimated CFT result x = | 
(dashed line). Right: The finite size linear fit (solid line) 
shows a mismatch with the estimated CFT result x ~ 4 . 



tonian for the spin-5 chain [1] consists of a bilinear and 
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where S a is a spin-S* operator acting on site a and 
periodic boundary conditions are imposed by sites 
parametrized as complex roots of unity rj a = exp (i^-a), 
a G 1,2, ... ,N. (Note that the biquadratic two-site term 
is contained in the three-site term in (1) as the special 
case f3 = 7.) For S = 1, this Hamiltonian has very 
recently been obtained independently through field theo- 
retical methods by Niclson, Cirac, and Sierra [24]. H1/2 
is the HSM [13, 14] (the three-site term trivially simpli- 
fies with the biquadratic term in this case). In Holstcin- 
PrimakofF representation, the ground state, which is the 
Gutzwiller wave function [25], takes the v = 1/2 Laugh- 

lin form ^ s (zi, . . . , z M ) = Yl^jizi-Zj) 2 Tjfli where 
M = vN and the z^s denote the coordinates of the spin 
flip operators acting on a spin polarized vacuum. The 
1st degree homogeneous polynomial multiplied with the 
squared Jastrow factor ensures that the state is both real 
and a spin singlet. 

<J>q s is the building block for all other spin-5 ground 
states. The spin-S 1 state is constructed by symmetriza- 
tion of the spin flip coordinates of k = 2S identical copies 
of X S^ S . Technically, this can be accomplished conve- 
niently in terms of auxiliary Schwinger boson creation 
operators and , where the symmetrization effectively 
reduces to a multiplication of k — 2S copies Vo) = 



(V>? s [at,&t]) I 0) [22]. The ground state energy yields 
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(1) is hence computed 
from the Schwinger projection of two HSM ground 
states and takes the explicit form ^Q faS (zi, 
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where M = N and 
and A is the 
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antisymmetrization operator. Note that as the HSM sin- 
glet ground state requires N even, the ground states for 
all S > 1 of (1) require N even as well. The case of N 
odd will be analyzed below. 

Conformal field theory — Let us now determine the 
long wavelength universality classes of "Hs- We employ 
the density matrix renormalization group (DMRG) [26] 
to compute the von Neumann entanglement entropy (EE) 
for large system sizes. To begin with, we checked the an- 
alytic formula for the ground state energy up to large 
values of S and N. For gapped models, the EE saturates 
to a constant value; for gapless models, the EE diverges 
logarithmically, and, for finite chains takes the analytic 
form [27] 
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FIG. 3. (color online) S = 1 energy spectra scaled to units 
Takhtajan-Babudjan Hamiltonian. Right: Hs for 5 = 1. 

The computation of the ground states in (1) turns out to 
be challenging. One reason are the long-range terms in 
the Hamiltonian which adds to the enlarged local basis 
for higher spin 5. For 5 = 1 (5 = 3/2) we have kept up 
to 800 (1200) states of the reduced density matrix and 
performed up to 12 (14) sweeps, which enabled us to keep 
the discarded entropy below 10 -7 (10~ 5 ). For (1) with 
5 = 1/2, 1, and 3/2, we confirm critical behavior due 
to the log-divergent behavior of the EE. For 5 = 1/2, 
1, and 3/2, we find central charges c — 1, c — 1.46(2) 
and c = 1.60(2) from chains up to L = 30 [Fig. 1], which 
are approximately consistent with the expected values 
cwzw = 3k/(k + 2) for WZWfc. The discrepancy be- 
tween the numerically obtained c and the asymptotically 
expected values is absent for 5 = 1/2 but sets in for 
higher 5. From finite size scaling, we cannot exclude 
the existence of non-monotonous corrections for small L 
spin-5 chains. Still, over all we find numerical indication 
that the low-energy theory of Hs is the SU(2) WZW25 
model within reasonable error bars. 

Logarithmic corrections — Having established that the 
CFT's related to the spin chains in (1) are of WZWfc 
type, we now address the question of logarithmic correc- 
tions. As stated before, the HSM as the 5 — 1/2 realiza- 
tion of (1) shows no finite size corrections as compared 
to the long wavelength limit. To analyze this issue in 
more detail, we calculate the scaling dimension x of the 
WZWfc primary fields, as this quantity is highly sensitive 
to finite size corrections [28]. Specifically, the value for 
the primary field at momentum tt can be extracted from 
E\ — Eq = 2ttvx/L, where v is the velocity parameter and 
E and Ei are the energies of the ground state and first 
excited state. Only small system sizes are needed to de- 
termine whether a spin chain resides at the conformally 
invariant fixed point or not. In Fig. 2 we have plotted 
Ei — Eq for %i/2 and Hi as computed by exact diago- 
nalization (ED) for chain lengths up to N = 16. (As Hi 
is not very sparse, the number of scattering elements is 
already of O(10 13 ) for N = 16.) With x = \ and § for 
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5 = 1/2 and 1 as predicted from CFT, we nicely observe 
the absence of finite size corrections for the HSM, but 
its presence for the S = 1 chain [Fig. 2]. This suggests 
that logarithmic corrections are present for 5 > g, and 
as a consequence that the models (1) do not represent the 
conformally invariant fixed points of WZWfc for k > 1. 

Spectrum of the S = 1 chain — We investigate the 
structure of the finite size spectrum of (1) for 5=1 
as obtained from ED. For this purpose, we compare 
it with the spectrum of the Takhtajan-Babudjan (TB) 
model [11, 12], which is likewise connected to WZW2 and 
describes the critical point between the dimerized phase 
and the Haldane gap phase of the bilinear-biquadratic 
spin-1 chain. The low energy spectra for N = 16 are 
plotted in Fig. 3. Aside from the singlet ground state at 
momentum k = 0, the spin multiplet quantum numbers 
of the lowest energy modes in the different momentum 
subsectors are related. The lowest level at momentum 
k = tt determines the finite size gap. Both spectra look 
very similar, showing a two-lobe feature [Fig. 3]. Still, 
the overlap of the Pfaffian and TB ground state is be- 
low .85 already for small system sizes, indicating that 
the models have different finite size structures. This be- 
comes explicit as we develop a specific trial state Ansatz 
for the 5=1 model in the following, which does not 
similarly apply to the TB model. 

Generalized Schwinger projection scheme — Up to now, 
we have only used our Schwinger projection scheme for 
HSM ground states to generate the exact N even ground 
states for the spin-5 models in (1). We now investigate 
what we can achieve when we employ the same approach 
for HSM excited states before projection. To begin with, 
we use the Schwinger projection scheme to obtain a suit- 
able trial state for the N odd ground state of the 5=1 
chain [Fig. 4]. Consider first the two 5=1/2 chains be- 
fore projection. For N odd, the low energy modes of the 
HSM are given by the single spinon branch which only 
covers one half of the momenta (—tt/2 < k sp < tt/2 and 
tt/2 < fc sp < 3tt/2 for N = 1 mod 4 and TV = 3 mod 4, 
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FIG. 4. (color online) (a) Construction of the S = 1 trial state Atrial (fc = 0) for N odd out of two Haldane-Shastry one-spinon 
states with momenta k = ±7r/2. (b) Spectrum of the S = 1 chain for JV = 15. The ground state at k — is ipo dd - ( c ) Overlap 
(Atrial | ^S dd ) vs. 1/JV. For TV ^ oo, the overlap is 0.9983(1). 



respectively) [Fig. 4a]. The ground state for N odd in 
the 5 = 1 chain is located at k = 0. (JV = 15 is shown in 
Fig. 4b.) As phases of the wave function are preserved 
in Schwinger boson notation, they multiply under pro- 
jection, and the total momentum of the projected wave 
function is given by the sum of single chain momenta be- 
fore projection. In order to find a good trial state , 4 7 Q dd 
for the 5=1 model (1), we choose the lowest energy 
one-spinon states in two HSM chains at — 7r/2 and 7r/2 
before projection. To fully match the quantum numbers 
of the target state, we also project the spin part onto the 
singlet component of the projected two one-spinon wave 
functions. With this construction, we obtain an excellent 
overlap with the N odd ground state of (1) with 5 = 1, 
which is of the order of 0.999 and hardly changes with 
system size [Fig. 4c] . Good overlaps can also be achieved 
for trial states to match other modes such as the lowest N 
even eigenstates of the 5=1 model in different momen- 
tum sectors [Fig. 3]. There, we project the lowest lying 
two-spinon HS eigenstate for different momenta and a 
HS ground state together. This correspondence explains 
how the lobe features for the 5=1 model are connected 
to the two-spinon levels of the HSM before projection. 

From there, a unified picture emerges: we can inter- 
pret the Schwinger boson projection of Haldane-Shastry 
states as the creation of spinon product states. This is 
a peculiar property of the Schwinger boson projection 
of Haldane-Shastry chains, as we found elsewhere that 
taking two copies of a generic 5=1/2 nearest neighbor 
Heisenberg model ground state as building blocks, for ex- 
ample, generates a trial ground state for a gapped spin-1 
model, where the spinons become confined [29]. In con- 
trast to this, the HSM projection still provides gapless 



states, which is somewhat intuitive as the HSM is a free 
spinon gas. As the Hilbert space after projection is signif- 
icantly smaller than the product space of the individual 
Haldane-Shastry chains before projection, this yields an 
overcomplete basis, and gives rise to selection rules that 
specify which many-spinons states before projection map 
onto each other after projection. This can be one way of 
defining the notion of a "blocking mechanism" connected 
to the non-Abelian statistics of the spinons [1]. Compar- 
ing the trial states we constructed here via Schwinger 
boson projection with the actual eigenstates of (1) for 
finite systems, we see that this construction yields rea- 
sonable approximations to the low energy modes of the 
system, even though it does not provide us with exact 
eigenstates. The trial state we have constructed for the 
5 = 1 model with N odd illustrates this point: since 
the spinons before projection have been chosen such that 
they cannot decay any further as they are located at the 
outer edges of the dispersion branches, we suppose that 
the "ideal" finite size WZW2 spin chain model in the 
sense of spinon product states should exactly correspond 
to this construction and give an overlap of unity. The ob- 
served deviation from that is another way to interpret the 
logarithmic corrections we found for the Hamiltonian (1) 
for 5 > 1/2. 
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